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The Laplacian gauge for gauge group SUiN) is discussed in perturbation theory. It is shown that 
to the lowest non-trivial order, ^^), configurations in the Laplacian gauge automatically satisfy the 
(finite difference) Landau gauge condition. Laplacian gauge fixed configurations are examined 
numerically and it is seen that to they do not remain in the Landau gauge. 



1. INTRODUCTION 

In the continuum, the Landau gauge is 
normally formulated as a differential condition. 
There are Gribov copies[l], of course, and these 
are seen on the lattice as well. On the lattice 
one normally formulates the Landau gauge 
condition as a maximization. 

Max ReTrY^ [U^(x) + ufix -fi)] 

g sites X ]i (1) 

U^(x) = g{x)U^{x)gHx^Vi) 

The reason for this is not so much to deal with 
Gribov copies as to eliminate a plethora of 
lattice artifact copies. These copies have the 
property that at each site the quantity to be 
maximized in Eq. (1) could be either a 
maximum, a minimum, or even a saddle. The 
global maximization condition defines a copy 
free gauge in principle, but there is no 
algorithm for finding such a global maximum. 
There are generically many local mtixima, 
which are the lattice versions of the continuum 
Gribov copies. 

In 1992, Vink and Wiese proposed a copy 
free gauge that could be implemented on the 
lattice, which they called the Laplacian 
gauge[2]. It shares with the Landau gauge the 
virtue of smoothness, so that it naturally 
relates gauge fields and other gauge variant 
quantities at a distance. Additional 
developments of this gauge have included a 
compact perturbative formulation of the 



Laplacian gauge for SU(2)[3], computations of 
the quark and gluon propagators[4], and 
studies of instantons[5]. 

In this tEilk, I will perturbatively apply Vink 
and Wiese's gauge fixing algorithm to generic 
configurations for an SU{N) or U(N) lattice 
gauge theory. Two principal results will be 
presented regarding the relationship between 
the Laplacian and the Landau gauges. One is 
that in the lowest non-trivial order, g^, 
configurations in the Laplacian gauge 
automatically satisfy the (finite difference) 
Landau gauge condition. The other is that 
ensembles of configurations numerically fixed 
to Laplacian gauge show a violation of the 
Landau gauge condition proportional to g^. 

2. THE LAPLACIAN GAUGE 

The Laplacian gauge is defined in terms of 
the lattice Laplacian (color indices suppressed) 

Mx,y) = Y [2 6,/ 

^ , (2) 

- 6 .U (x) 6 . uHy)] 

A is a positive Hermitean operator. For SUiN), 
one forms an NxN matrix field from the 
eigenvectors with the N lowest eigenvalues, 

M/x) = (fSx))j (3) 

performs a polar decomposition at each site, 
and defines a gauge transformation Q{x) as the 
adjoint of the unitary part of M(x). Applying 



this gauge transformation to the lattice hnk 
variables transforms the lattice to the 
Laplacian gauge. Vink and Wiese showed that 
the Laplacian gauge was free of Gribov copies, 
that is, the transformation to Laplacian gauge 
gave the same result, wherever one started on 
the gauge orbit of a given configuration. This 
does not violate Singer's theorem[6] that no 
continuous (in the field variables) gauge could 
be unique. There are discontinuities in the 
Laplacian gauge; they occur at configurations 
where A(x^) has degeneracies among its first N 
eigenvalues, which signals level crossing. 

3. LATTICE PERTURBATION THEORY 

Let us expand the Laplacian gauge 
algorithm in a power series in the bare lattice 
coupling. Define a lattice gauge potential by 

U^ix) = (4) 

and use the power series for to expand the 
Laplacian 

Mx,y) = A<o>(x,y) 

^gA^'Kx,y) ^g^A^^Kx,y) + • • 

In zeroth order, i.e., U^ix) =1^, the Laplacian 
has an iV-fold degenerate zero eigenvalue, 
which, since A"'' is also a positive operator, are 
the lowest eigenvalues. Their eigenfunctions 
are the N linearly independent constant 
AT-vectors, which we denote as 

= ?f i = l,2,...,N (6) 

We call the matrix formed from them H"'' . 

From the orthogonality of the /"/°'(x) and 
their constancy in x we can choose their 
normalizations so that H"^' is unitary. The 
degeneracy makes it easy to develop the 
analysis, but it also means that the 
perturbative expansion of the Laplacian gauge 
is a case of degenerate perturbation theory, in 
which higher order effects determine which 
linear combinations of the functions are the 
proper zeroth order eigenfunctions, i.e. remain 
unmixed in higher order. This higher order 
information, which in fact appears in second 
order, is needed even to correctly construct the 



zeroth order Laplacian gauge transformation^ 
Through second order, the equations for the 
eigenfunctions are 

^(0) ^(1) + ^(1) ^(0) ^ ■^m ^(0) 

(7) 

A<°' /■.*^* + A<" + A<2' fl^^ = A." /■/" + X^^^ fl°^ 

From the first equation it is seen that the first 
order shift in the eigenvalues vanishes, and the 
first order wave functions are a simple linear 
transformation of the zeroth order ones. This 
gives for the first order matrix M 

M^^Kx) = -iA*"' 'u,y)A^A^(3/)H<o> 
where 

AA(^) = E[A,(y)-V3'-A)] (9) 

is the lattice divergence of the gauge potential 
based on Eq.(4), and A**^^ is the inverse of the 
free Laplacian restricted to the space 
orthogonal to the f-°\x) . 

Projecting the second order equation on the 
space of constant functions gives an ordinary 
NxN matrix eigenvalue problem: 

p(0)(^(2) _ a<i>a"''"^A<i>)/'/'^' = Xffl'^^ (10) 

Its solutions, which break the degeneracy 
between the N lowest eigenvalues, are the 
stable zeroth order solutions. The orthogonal 
projection of this equation gives the second 
order term in M: 

M<2)(x) = -A*°*"'(a*2' - A*"A*"*"'a*")h<«> 

-Z2(x)H<«> ^ ^ 

These operator exercises are easily carried out 
in momenttim space. The further construction 
of the transformation to Laplacian gauge 
proceeds one site at a time in coordinate space. 
It is straightforward to check that Z^(x) is 

Hn the case of SU(2), the degeneracy persists to all orders 
because of conjugation symmetry, i.e., because the 
fundamental representation of SU(2) is unitarily equivalent 
to its complex conjugate. This allows one to choose the 
matrix M(x) proportional to the unit matrix[3]. 



antihermitian, so that to first order M(x) is a 
unitary matrix at each site. However, at 
second order its unitary part must be exphcitly 
extracted, giving through second order 

Qt(x) = Mix) \M\x)M{x)Y^'^ 



w(0) (12) 



3.1 Transformation of the gauge potential 

It is worthwhile to express the 
transformation to Laplacian gauge in terms of 
the gauge potential, using the expansion of 



-igA (x) 



Q(x)[7^(x)Q^(x + iii) (13) 



generated with the simple Wilson action for the 
gauge group SU(3) at values of (3=6/^^ 
corresponding to g = .1, .2, ... ,.9, 1 . After fixing 
to Laplacian gauge, we evaluated the lattice 
divergence of the standard lattice gauge 
potential. 



A^ix) 



[U(x) - U!ix)] 



Traceless 



2ig 



(19) 



The average of over each ensemble is 

plotted vs the bare lattice coupling in Fig. 1. 

Fig. 1 Deviation of Laplacian Gauge 
Fixed Lattices From Landau Gauge 



Because of the interplay between factors of^ in 
the expansions of Q and U^, the coefficient of 
g^, which is the first non-trivial term, is the 
term in the gauge potential transformed 
to Laplacian gauge. Explicitly, it is 

A^^%) = H«»^[A^(x) - il.ix) + t;,(x +A)]S<°'(14) 

As has been observed from the first discussions 
of the Laplacian gauge, the lattice divergence of 
this gauge potential vanishes 



-(0)t 



7(0) 



(15) 







The new result presented here, eliding much 
algebra, is that the next term in the expansion 
of the gauge-fixed gauge potential, 

. _ 1 ,-.(( 

~~ 2 



-i(Z2-Zj) + i(Z2<^^*-Z<^^*^)}H(0' 



(16) 



where 



(17) 



also satisfies the Landau gauge condition 



4. NON-PERTURBATIVE FIXING 



(18) 



The foregoing considerations can be verified 
in simulations. To do this, we applied the 
Laplacian gauge fixing procedure to ten 
ensembles of 20 well thermalized 6* lattices 
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The curve is a quadratic fit to the square root of 
the ensemble average of TriA^^)^ averaged 
over lattice sites. 

The graph shows two relevant results. One 
is that the Landau gauge condition clearly fails 
to second order in perturbation theory. The 
other is that even at ^ = 1, which corresponds 
to p = 6, the magnitude of A^/l^ is only about 
1% of its value on unfixed lattices. 

REFERENCES 

1. V.N. Gribov, Nucl. Phys. B 139 (1978) 1 

2. J.C. Vink and U.^. Wiese, Phys. Lett. B 289 (1992) 122; 
J. C. Vink, Phys. Rev. D 51 (1995) 1292 

3. P. van Baal, Nucl. Phys. B (PS) 42 (1995) 843 

4. C. Alexandrou, Ph. de Forcrand, and F. FoUana, hep- 
lat/0008012, hep-lat/0009003; O. Philipsen, hep- 
lat/0106006; P.O. Bowman, U.M. Heller, and A.G. 
Williams, hep-lat/0 110081 

5. Ph. de Forcrand and M. Pepe, Nucl. Phys. B (PS) 94 
(2001) 498 

6. I.M. Singer, Comm. Math. Phys. 60 (1978) 7 



